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The preliminary CERN SPS NA49 Pb+Pb 158 GeV/A one- and two-particle /i~-spectra 
at mid-rapidity are consistent with a source of temperature T « 130 MeV, lifetime r ~ 9 
fm/c, transverse flow rjf ~ 0.35, and a transverse geometric size which is twice as large 
as the cold Pb nucleus. 
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00 ■ 1. Introduction 

Hadronic one- and two-particle spectra provide for each particle species information 
about the phase-space distribution S(x,K) of the hadronic emission region [l|]. Once 
reconstructed from the measured spectra, S(x, K) allows to distinguish between different 
dynamical scenarios of heavy ion collisions. It provides an experiment-based starting point 
for a dynamical back extrapolation into the hot and dense early stage of the collision, 
where quarks and gluons are expected to be the relevant physical degrees of freedom. 
Our reconstruction of S(x, K) is based on the hadronic one- and two-particle spectra 
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The triple-differential hadronic one-particle spectrum (P tests the momentum-dependence 
of S(x, K) only. Its azimuthal 0-dependence with respect to the reaction plane is 
parametrized by the harmonic coefficients v n [0 . Information about the space-time struc- 
ture of S(x, K) can be obtained from the relative momentum dependence of the two- 
particle correlator C(K, q), q = Pi — P2- This q-dependence is usually parametrized 
via the Hanbury-Brown Twiss (HBT) radius parameters _R^(K) which depend on the 
average pair momentum K = |(pi + p2)- Depending on the Gaussian parametrization 
adopted, the indices i, j 'in (Q) run either over the Cartesian directions long (parallel to the 
beam), out (parallel to the transverse component Kj_) and side, or over the corresponding 

Yano-Koonin coordinates q± = w ql + gf, ^3 and q° = q • K/_K" 0- 
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Figure 1. LHS: \ 2 contour plot of a fit to the NA49 P /i _ -spectrum. Dashed lines are for 
constant values of tA/T. RHS: different combinations of temperature T and transverse 
flow rjf account for the same one-particle slope. 



2. Reconstructing an azimuthally symmetric source 

A typical data analysis starts from a simple ansatz for the phase space distribution 
S(x, K) in terms of very few, physically intuitive fit parameters, 



S„(x,p) = S^ T (x,p) +J2 S R^( X ^P) ' ( 3 ) 

R 

This model e.g. assumes local thermalization at freeze-out with temperature T within 
a space-time region of transverse Gaussian width R, emission duration At, longitudinal 
extension tqAt], where 77 = hki[(t + z)/(t — z)], and average emission time r . The 
model allows for dynamical source correlations via the hydrodynamic flow 4- velocity u^{x). 
We assume a linear transverse flow profile r] t (r) = Vf with variable strength rjf, 
and Bjorken scaling of the flow component in the longitudinal direction, vi = z/t, rji = 
I ln[(l + v{)/ (1 — vi)] = 1]. Resonances are produced in thermal abundances with proper 
spin degeneracy 2 Jj + 1 for each particle species i. Their contribution to the pion yield is 
obtained by propagating them along their classical path = + according to an 
exponential decay law ||, 

S R ^(x,p) = J R Jd*X JdrTe- rT 5^[x-(x + { I r)]Sf(X,P), (5) 

where / R is the integral over the available resonance phase space for isotropic decays. We 
include all pion decay channels of p, A, K*, £*, u, rj, 77', Kg, S and A with branching 
ratios larger than 5 percent. The model parameters T, r]f, R, Ar), At, tq can then be 
determined via the following strategy: 

1. transverse one-particle spectrum dN/dM\ determines blue-shifted temperature T s q: 
The slope of dN/dM]_ is essentially given by T e fj = TJ(1 +r?/)/(l —77/). Hence, 



different combinations of T and rjf can account for the same data, see Fig. [1]. 

dN/dM]_ and R±(M±) disentangle temperature and transverse flow. R± fixes trans- 
verse extension R: 
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Figure 2. Yano-Koonin-Podgoretskii HBT-radius parameters. The preliminary NA49 
Pb+Pb data || are for h + h + (squares) and h~h~~ (diamonds) correlations. 



The M ± -slope of R ± is proportional to nj/T, R 2 ± & R 2 /(l + M ± r]j/T) §,§. This 
slope allows in combination with dN/dM\ to specify T and 7]f, see Fig. |2|. The 
overall size of R± determines the Gaussian width R. We find 77/ ~ 0.35, T « 130 
MeV, R « 7 fm. 

3. Ru determines r Q , width of rapidity distribution dN/dY determines An: 

In principle, R\\ depends on tq, An and At 0. We have fixed An = 1.3 by matching 
the width of the pion rapidity distribution. The data then favour clearly a lifetime 
of To ~ 9 fm/c (see Fig. |2|) but are not very sensitive to the emission duration At ||. 
The present plots are obtained with At = 1.5 fm/c. 

4. R discards opaque sources: 

For the model (|3|-[5]), the YKP-parameter Rq is mainly sensitive to the temporal 
aspects of the source. The large statistical uncertainties for i? do not allow to 
constrain the model parameter space further, see Fig. ||. Models of opaque sources 
including an opacity factor in (|) are excluded already by the present data ||. 

The radius of a cold Pb nucleus corresponds to i?pb(cold) ~ 3.5 fm in the Gaussian 
parametrization (f|), i.e., the experimental data indicate a very large source Rm2 _R Pb (cold) 
This is dynamically consistent with a collision scenario in which the initially produced 
pressure gradients result in a significant transverse flow r]f = 0.35, driving the expansion 
of the system over a time of 9 fm/c to twice its initial size. The temperature may have 
decreased substantially during this expansion; the data indicate 130 MeV at freeze-out. 
These conclusions are further supported by the analysis presented by G. Roland |§, which 
is based on approximate analytical formulas. We have extracted the model parameters by 
comparison to numerical model calculations, following the above strategy. They are not 
obtained from a simultaneous fit to all observables, and hence we do not quote errors. 

3. Particle interferometry for collisions with finite impact parameter 

The reaction plane analysis of triple-differential one-particle spectra ([!]) has been dis- 
cussed extensively in this conference e.g. by A. Poskanzer, J.-Y. Ollitrault, and S. 
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Voloshin. The general strategy for linking this analysis to an azimuthally sensitive parti- 
cle interferometry is based on the ^-dependence of the HBT-radii |T(J , where $ measures 
the azimuthal angle of Ki relative to the reaction plane, 

oo oo 

Rl(K ± , $, Y) = R 2 j)0 (K ± , Y)+2 £ R% 2 [K U Y) cosn$+2 £ Rt j>n 2 (K ± , Y) sinn$ .(6) 

n=l n=l 

Here, the many harmonic coefficients i?*- # make a direct comparison to experimental data 
impossible. However, various relations hold amongst these coefficients, since the leading 
anisotropy in realistic source models S(x, K) can be quantified by very few parameters 
only. Using the symmetries of the system and assuming that elliptic deformations domi- 
nate we find |H| 

D srf R c 2 eh R s 2 srf f? c 2 ssi R s 2 srf R c 2 srf f? s 2 m > 1 

w "i,m J W,m ±L ol,m ±L ol,m ±L sl,m ±L sl,m i " L — x i \' J 

«i « ^,i 2 « ^,i 2 « > (8) 

«2 « -^, 2 2 « ^/ « • (9) 

A violation of Eqs. (B|)-(|5|) by experiment would indicate strong higher order deformations 
and rule out many model scenarios. On the basis of Eqs. (0)-©! an azimuthally sensitive 
parametrization of the two-particle correlator involves only two additional fit parameters, 

C^(K, q) « 1 + A exp[-i? Oj0 2 q 2 - R s / q 2 - Ri, 2 qf - 2 R i fi 2 q qi] 

x exp [ - at (3 q 2 Q + q 2 ) cos($ - ijj R ) + 2a x q Q q s sin($ - ijj R )\ 
x exp [ - a 2 (q 2 - q 2 ) cos 2($ - tp R ) + 2 a 2 q q s sin 2($ - VflKJlD) 

The anisotropy parameter a,\ vanishes at mid-rapidity or if the source contains no dynam- 
ical correlations. It characterizes anisotropic dynamics. The parameter a 2 characterizes 
the elliptic geometry. The parameters a.\ and a 2 can be determined from event samples in 
spite of the uncertainty in the eventwise reconstruction of the reaction plane. For details, 



see Ref. p] 
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